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ABSTRACT 

In order to test the reliability of determinations of the shapes of galaxies' dark 
matter halos, we have made such measurements for the Milky Way by two independent 
methods. First, we have combined the measurements of the over-all mass distribution 
of the Milky Way derived from its rotation curve and the measurements of the amount 
of dark matter in the solar neighborhood obtained from stellar kinematics to determine 
the flattening of the dark halo. Second, we have used the established technique based 
on the variation in thickness of the Milky Way's Hi layer with radius: by assuming 
that the Hi gas is in hydrostatic equilibrium in the gravitational potential of a galaxy, 
one can use the observe flaring of the gas layer to determine the shape of the dark 
halo. 

These techniques are found to produce a consistent estimate for the flattening 
of the dark matter halo, with a shortest-to-longest axis ratio of q ~ 0.8, but only 
if one adopts somewhat non-standard values for the distance to the Galactic centre, 
i?o, and the local Galactic rotation speed, Qq- For consistency, one requires values 
of i? Ss 7.6 kpc and Oo ^ 190 km s -1 . The results depend on the Galactic constants 
because the adopted values affect both distance measurements within the Milky Way 
and the shape of the rotation curve, which, in turn, alter the inferred halo shape. 
Although differing significantly from the current IAU-sanctioned values, these upper 
limits are consistent with all existing observational constraints. 

If future measurements confirm these lower values for the Galactic constants, then 
the validity of the gas layer flaring method will be confirmed. Further, dark matter 
candidates such as cold molecular gas and massive decaying neutrinos, which predict 
very flat dark halos with q ^ 0.2, will be ruled out. Conversely, if the Galactic constants 
were found to be close to the more conventional values, then there would have to be 
some systematic error in the methods for measuring dark halo shapes, so the existing 
modeling techniques would have to be viewed with some scepticism. 

Key words: Galaxy: structure - Galaxy: kinematics and dynamics - Galaxy: solar 
neighbourhood - Galaxy: fundamental parameters - Galaxy: stellar content - ISM: 
general 



1 INTRODUCTION 

The speed at which spiral g alaxies rotate remains relatively 
const ar 



I otri lu laigu tadii (e.g., Rubin, Fbid fe ThuiiiiaiiJ 



dieted by gal axy formation models in which hot dark matter 
is dominant (Peebles 1993). Cosmological cold dark matter 
simulations typically result in triaxial dark halos ( [Warren 



198C; B oama 1981), which implied that they contain large 



et al. 1992) while the inclusion of gas dynamics in the sim- 



amounfc of dark mat Ler. The nature of Lilly malerial remalny 

obscure 



ulations resul t s in somewhat flattene d , oblate halos ( Katz 



the dark matter halo, as quantified by its shortest-to-longest 
axis ratjlo, q = u/a. The roundest lialoy with q — 0.8 arc pre 



Cnnu l')')lj; Judiff fe Martinet 1994| ; [Dubinski 1994| ) with 
candidates such a s cold molecular gas (Pfenniger, Cornbee 



buL one key diagnostic is pi uvided by the shape uf - . 5 ±0.15 ( fcubinski 1994| ). Models' fo r other dark matter 
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fe Martinet 1994 ) and massive decaying neutrinos ( jciama 
199C) require halos as flat as q ~ 0.2. Clearly, the determi- 
nation of q for real galaxies offers a valuable test for discrim- 
inating between these cosmological models. 
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Figure 1. A summary of the existing estimates of galaxies' dark 
matter halo shapes as parameterized by their shortest-to-longest 
axis ratio, q. For each individual galaxy we indicate the identifi- 
cation and techni que used to determine the halo's flattening: 1) 
Flar ing g as layer (Oiling 1996b|_ Bccquacrt. Combes &, Viallefond 



NGC 4650A 



A0136-0801 



AM2020-504 



f997); 1) Warping gas layer (|Hofn er & Sparke 1994); 3) X-ray 



isophot es, ( Buote fc Canizares 199£ ): 4 j Polar ring galaxies ( Arn 
aboldi c \ al. 1993 ; |Sackctt et al. 199-1 :_ Sackctt. fe 1'oggc 199' ); 5) 



prec essin g dusty disk (Stciman-Cameron, Kormendy Sz Duriscn 
1992 [). ^ he dotted line shows the predicted distribution for the 
shapes of halos in a cold dark matter cosmology (q = 0.5 ± 0.15; 
Dubinski 1994). 



In Fig. [l], we summarize the existing estimates of halo 
shape, derived using a variety of techniques. It is evident 
from this figure that there is only a very limited amount of 
data available for measuring the distribution of q in galax- 
ies. Rather more worrying, though, is the fact that different 
techniques seem to yield systematically different answers. 
The warping gas layer method, in which the shape of the 
halo is inferred by treating any warp in a galaxy's gas layer 
as a ben ding mode in a flattened potential (rlofner & Sparke 
199 4]), | eems to imply that dark halos are close to spheri- 
cal, with q <; 0.8. Conversely, the flaring gas layer tech- 
nique, which determines the halo shape by assuming that 
the gas layer is in hydrostatic equilibrium in the galaxy's 
gravitational potential and uses the thickness of the layer 
as a diagnost ic for the distribution of mass in the galaxy 
(Oiling 1996b), yields much flatter halo shape estimates with 
q ^ 0.4. Although the numbers involved are rather small, 
there do seem to be real differences between the results ob- 
tained by the different methods. Thus, either these tech- 
niques are being applied to systematically different classes 
of galaxy, or at least some of the methods are returning 
erroneous results. 

In order to determine which techniques are credible, we 
really need to apply several methods to a single galaxy, to 
see which produce consistent answers. As a first step towards 



such cross-calibration, this paper compares the shape of our 
own galaxy's dark matter halo as inferred by two distinct 
techniques: 

(i) Stellar kinematics. Our position within the Milky Way 
means that we have access to information for this galaxy 
that cannot be obtained from other systems. In particular, 
it is possible to measure the total column density of material 
near the Sun using stellar kinematics. After subtracting off 
the other components, we can infer the local column density 
of dark matter. By comparing this density close to the plane 
of the Galaxy to the over-all mass as derived from its rate 
of rotation, we can obtain a direct measure of the shape of 
the halo. 

(ii) The flaring gas layer method. As outlined above, this 
technique assumes that the Hi emission in the Milky Way 
comes from gas in hydrostatic equilibrium in the Galactic 
potential, from which the shape of the dark halo is inferred. 
Since the results of previous applications of this method give 
results somewhat out of line with the other techniques, it is 
important to assess the method's credibility. 

As well as providing a check on the validity of the flaring 
gas layer method, these analyses will also provide another 
useful datum on the rather sparsely populated Fig. |l| 

The remainder of the paper is laid out as follows. Both 
of the above methods rely on knowledge of the Milky Way's 
rotation velocity as a function of Galactic radius - its "rota- 
tion curve" - so Section ^ summarizes the data available for 
estimating this quantity, and the dependence of the inferred 
rotation curve on the assumed distance to the Galactic cen- 
tre and local rotation velocity. The analysis of the shape of 
the dark halo requires that we decompose the Milky Way 
into its visible and dark matter components, so Section ^ 
discusses the construction of a set of models consistent with 
both the photometric properties of the Galaxy, and its mass 
properties as inferred from the rotation curve. In Section ^ 
we show how these mass models can be combined with the 
local stellar kinematic measurements to determine the shape 
of the dark halo. Section ^presents the application of the gas 
layer flaring technique to the mass models. Section ^ com- 
bines the results derived by the two techniques and assesses 
their consistency. The broader conclusions of this work are 
drawn in Section ^ 



2 THE OBSERVED ROTATION CURVE 

Our position within the Milky Way complicates the geome- 
try when studying its structure and kinematics. It is there- 
fore significantly harder to determine our own galaxy's ro- 
tation curve, Q(R), than it is to derive those for external 
systems. More directly accessible to observation than Q(R) 
is the related quantity 



W(R) = Ro 



R Ro 



(1) 



where Rq and Oo are the distance to the Galactic centre 
and the local circular speed. If one assumes that material 
in the galaxy is in purely circular motion, then some simple 
geometry shows that, for an object at Galactic coordinates 
{l,b} with a line-of-sight velocity v\ os , 
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W 



I-'lo 



sin I cos b 



(2) 



(Binney & Merrifield 1998, §9.2.3). Thus, if one measures 
the line-of-sight velocities for a series of objects at some 
radius R in the Galaxy, one has an estimate for W(R). By 
adopting values for Ro and ©o, one can then use equation (M 
to determine the rotation speed at that radius. 

In practice, the difficulty lies in knowing the Galactic 
radii of the objects one is looking at. One solution is to look 
at standard candles, whose distances can be estimated, and 
hence whose radii in the Galaxy can be geometrically de- 
rived. Alternatively, one can select the subset of some tracer 
- usually Hi or H2 gas - whose line-of-sight velocities and 
Galactic coordinates places it at the same value of W , and 
hence at the same radius. One can then use the proper- 
ties of this cylindrical slice through the Galaxy to infer its 
radius. For example, in the inner Galaxy, all the material 
in a ring of radius R will lie at Galactic longitudes of less 
than i max = sin _1 (i?/i?o), so one can use the extent of the 
emission on the sky of each Vy-slice to infer its radius - an 
approach traditionally termed the "tangent point method" 
[see, for example, Malhotra (1994,1995)]. At radii greater 
than Ro this method is no longer applicable as the emission 
will be visible at all values of I. However, by assuming that 
the thickness of the gas layer does not vary with azimuth, 
one can use the observed variation in the angular thickness 
of the layer with Galactic longitude to estimate the radii of 
such W slices (Merrifield 1992). 

For the remainder of this paper, we use Malhotra's 
(1994, 1995) tangent point analysis to estimate W(R) in 
the inner Galaxy. For the outer Galaxy, we have combined 
Merrifield's (1992) data with Brand & Blitz's (1993) anal- 
ysis of the distances to HII regions, from which standard 
candle analysis the rotation curve can be derived. 

In order to convert W(R) into O(-R) using equation (|l|), 
we must adopt values for the Galactic constants, Ro and 
Oo. Unfortunately, there are still significant uncertainties 
in these basic parameters. In the case of Ro, for example, 
the extensive review by Reid (1993) discussed measurements 
varying between Ro — 6.9 ± 0.6 kpc and Ro = 8.4 ± 0.4 kpc. 
Even more recently there have been few signs of conver- 
gence: Layden et al. (1996) used RR Lyrae stars as stan- 
dard candles to derive Ro = 7.2 ± 0.7 kpc, while Feast 
& Whitelock (1997) used a Cepheid calibration to obtains 
i?o = 8.5 ± 0.5 kpc. The constraints on Oo are similarly 
weak: a recent review by Sackett (1997) concluded that a 
value somewhere in the range Oo = 210 ± 25 kms -1 pro- 
vided the best current estimate. It should also be borne in 
mind that the best estimates for Ro and Oo are not in- 
dependent. Analysis of the local stellar kinematics via the 
Oort constants gives quite a well-constrained value for the 
ratio O /-Ro = 26.4 ± 1.9 km s" 1 kpc" 1 (Kerr & Lynden- 
Bell 1986), so a lower-than-average value of Ro is likely to 
be accompanied by a lower-than-average value for Oo. Cur- 
rently, the best available measure of the local angular ve- 
locity of the Milky Way is based on VLBI proper motion 
measurements of SgrA*. Assuming that SgrA* is at rest 
with respect to the Galactic centre, Reid et al. (1999) find 
©0/-R0 = 27.25 ± 2.5km s _1 kpc -1 , consistent with the value 
proposed by Kerr & Lynden-Bell. 

To illustrate the effect of the adopted values of Oo and 
Ro on the derived rotation curve, Fig. ti shows O(-R) for 



300 



200 



100 



300 



200 



100 




5 10 15 20 

Galactocentric Radius [ kpc ] 



Figure 2. The rotation curve for the Milky Way for val- 
ues of Ro = 7.1kpc, ©0 = ISSkms" 1 , and Ro = 8.5kpc, 
©0 = 220 kms" 1 . The figure also shows one of the ways in which 
the rotation curve can be decomposed into the contributions from 
different mass components: the bulge (dotted line); the stellar disk 
(filled circles); the Hi layer (crosses, where negative values mean 
that the force is directed outwards); the H2 layer (circles); and 
the dark halo (dashed line). The best fit model, which is obtained 
by summing the individual components in quadrature, is shown 
as a full line. 



two of the more extreme plausible sets of Galactic parame- 
ters. Clearly, the choice of values for these quantities affects 
such basic issues as whether the rotation curve is rising or 
falling in the outer Galaxy. Given the current uncertainties, 
it makes little sense to pick fixed values for the Galactic 
constants, so in the following analysis we consider mod- 
els across a broad range of values, 5.5 kpc < Ro < 9 kpc, 
165 kms" 1 < Oo < 235 kms -1 . 



3 MASS MODELS 

In order to relate the rotation curve to the shape of the 
dark halo, we must break the gravitational potential respon- 
sible for the observed 0(i?) into the contributions from the 
different mas s components. As is usually don e in this de- 
composition (e.g. van Albada fc Sancisi 1981: ; Kent 1987 ; 
Bcgcman 1989]; |Lake fc Feinswog 1989t |Broeils 1992| ; piling 



1995 



Oiling 1996c; Dehnen & Binney 199c), we adopt a 



model consisting of a set of basic components: 
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(i) A stellar bulge. Following Kent's (1992) analysis of 
the Galaxy's K-band light distribution, we model the Milky 
Way's bulge by a "boxy" density distribution, 



pb{R, z) <x Ko(s/hb), where s 4 = R 4 + (z/qb) 4 



(3) 



This modified Bessel function produces a bulge that appears 
exponential in projection. The observed flattening of the K- 
band light yields gj, = 0.61, and its characteristic scalelength 
is hb = 670 pc (Kent 1992). The constant of proportionality 
depends on the bulge mass-to-light ratio, Tf,, which we leave 
as a free parameter. 

(ii) A stellar disk. The disk is modelled by a density dis- 
tribution, 



p d (R,z) <x exp(-R/h d )sech(z/2z d ). 



(4) 



The first term gives the customary radially-exponential disk. 
The appropriate value for the scalelength, hd, is still some- 
what uncertain - Kent, Dame & Fazio (1991) estimated 
hd — 3 ± 0.5 kpc, while Freudenreich (1998) found a value 
of 2.5 kpc. We therefore leave this parameter free to vary 
within the range 2 kpc < hd < 3 kpc. The ^-dependence 
adopts van der Kruit's (1988) compromise between a sech 2 
isothermal sheet and a pure exponential. For simplicity, we 
fix the scaleheight at Zd = 300 pc. However, the exact in- 
dependence of pd was found to have no impact on any of 
the following analysis. Once again, the constant of propor- 
tionality depends on the mass-to-light ratio of the disk, Td, 
which we leave as a free parameter. 

(iii) A gas disk. From the Hi data given by Burton (1988) 
and Malhotra (1995) and the H2 column densities from 
Bronfman et al. (1988) and Wouterloot et al. (1990), we 
have inferred the density of gas as a function of radius in 
the Galaxy. This analysis treats the gas as an axisymmet- 
ric distribution, and neglects the contribution from ionized 
phases of the interstellar medium, but does include a 24% 
contribution by mass from helium (Olive & Steigman 1995) . 

(iv) A dark matter halo. We model the dark halo as a 
flattened non-singular isothermal sphere, which has a den- 
sity distribution 



p h (R,z) = p h - 



Rh 



(5) 



'Rl+R^ + iz/q) 2 ' 

where pu is the central density, Rh is the halo core radius, 
and q is the halo flattening, which is the key parameter in 
this paper. 

The procedure for calculating a mass model from these 
components is quite straightforward. For each pair of Galac- 
tic constants, Ro and Oo, we vary the unknown parameters 
{Tb, hd, Yd, ph, Rh, q} to produce the mass model that has 
a gravitational potential, <fr(-R, 2), such that 



m - ( *1 



1/2 



(6) 



provides the best fit (in a minimum \ 2 sense) to the observed 
rotation curve, Q(R). Examples of two such best-fit models 



As is well known (3. 


g. van 


Albada & Sancisi 1986 


Kent 1987; Begeman 1989 


Lake 


& Feinswog 198S; Broeils 


1992; Oiling 1995; Oiling 


1996c; 


Dehnen & Binney 1998), 



so many different combinations of the individual compo- 
nents can reproduce the observed rotation curve with almost 
equal quality of fit. We have therefore searched the entire 
{Tb, h d ,T d, ph,Rh, q} parameter space to find the complete 
subset of values that produce fits in which the derived value 
of x 2 exceeded the minimum value by less than unity. 



4 HALO FLATTENING FROM LOCAL 
STELLAR KINEMATICS 

Although the analysis of the previous section tells us some- 
thing about the possible range of mass models for the Milky 
Way, it does not place any useful constraint on the shape of 
the dark halo: for any of the adopted values of Ro and Oo, 
one can find acceptable mass models with highly-flattened 
dark halos (q ~ 0), round dark halos (q ~ 1), and even 
prolate dark halos (q > 1). We therefore need some further 
factor to discriminate between these models. 

One such constraint is provided by stellar kinematics in 
the solar neighbourhood. Studies of the motions of stars in 
the Galactic disk near the Sun imply that the total amount 
of mass within 1.1 kpc of the Galactic plane is S1.1 = (71 ± 
6) Mrs, PC" 2 (Kuijken & Gilmore 1991). Clearly, the value of 
S1.1 provides an important clue to the shape of the Milky 
Way's dark halo; in general, a model with a highly-flattened 
dark halo will place a lot of mass near the Galactic plane 
leading to a high predicted value for £1.1, while a round 
halo will distribute more of the dark matter further from 
the plane, depressing £1.1. 

The dark halo is not the only contributor to £1.1. In 
particular, the stellar disk has a surface density in the solar 
neighbourhood, £» , which may contribute a significant frac- 
tion of £1.1. There must therefore be a fairly simple relation 
between the adopted value of £* and the inferred halo flat- 
tening, q. Specifically, as one considers larger possible values 
of £», the amount of dark matter near the plane must de- 
crease so as to preserve the observed value of £1.1. Such a 
decrease can be achieved by increasing q, thus making the 
dark halo rounder. 

This inter-relation is illustrated in Fig. ^[ Here, for one 
of the sets of possible Galactic constants, we have considered 
all the mass models that produce an acceptable fit to the 
rotation curve and reproduce the meaured value of £1.1. 
For each acceptable model, we have extracted the value of 
the halo flattening q and the mass of the stellar disk in the 
solar neighbourhood, £*. For the reasons described above, 
these quantities are tightly correlated, and Fig. ^ shows the 
linear regression between the two. 

Clearly, we have not yet derived a unique value for the 
halo flattening: by selecting models with different values of 
£* , we can still tune q to essentially any value we want. How- 
ever, £* is not an entirely free parameter. From star-count 
analysis, it is possible to perform a stellar mass census in 
the solar neighbourhood, and thus determine the local stel- 
lar column density. Until relatively recently, this analysis has 
been subject to significant uncertainties, with estimates as 
high as £* = 145M Q pc~ 2 (|Bahcall 1984]). However, there 



such mass decompositions are by no means unique: there is 
near degeneracy between the various unknown parameters, 



is now reasonable agreement between the various analyses, 
with more recent published values of 35 ± 5Mq pc -2 (Kui- 
jken & Gilmore 1989) and 26 ± 4M© pc" 2 (Gould, Bahcall 
& Flynn 1997). If we adopt Kuijken & Gilmore's slightly 
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Figure 3. The flattening of the model Galactic halo, q, as a 
function of the adopted model's column of stellar disk mass in the 
Solar neighbourhood, E* . We show the best fits to the stellar kine- 
matic constraints (dotted line) and the flaring of the Hi layer (full 
line). The shaded region shows values of S* that lie more than 
1-CT from the observed value. This particular set of calculations 
has been made assuming Ro = 7.1 kpc and 0o = 185kms — 1 . 



more conservative error bounds on E», the shaded regions 
in Fig. ^| no longer represent acceptable models as they pre- 
dict the wrong local disk mass density, so we end up with a 
moderately well-constrained estimate for the halo flattening 
of q = 0.7 ±0.1. 

Note, though, that although this analysis returns a good 
estimate for q, Fig. ^ only shows the value obtained for a par- 
ticular set of Galactic constants. As discussed above, chang- 
ing the adopted values for Ro and Oo alters the range of 
acceptable mass models, which, in turn, will alter the de- 
rived correlation between E* and q. The uncertainties in the 
Galactic constants are still sufficiently large that the abso- 
lute constraints on q remain weak. Nevertheless, it is to be 
hoped that the measurements of Ro and Go will continue to 
improve over time, leading to a unique determination of q 
by this method. Further, as we shall see below, it may be 
possible to attack the problem from the other direction by 
using other estimates of q to help determine the values of 
the Galactic constants. 



from its rotation curve, while the degree to which this mass 
distribution is flattened is derived using the properties of a 
tracer population close to the Galactic plane. In this case, 
the tracer is provided by the Hi emission from the Galactic 
gas layer. The thickness of the Milky Way's gas layer is dic- 
tated by the hydrostatic balance between the pull of gravity 
toward the Galactic plane and the pressure forces acting 
on the gas. As the density of material in the Galaxy de- 
creases with radius, the gravitational force toward the plane 
becomes weaker, so the equilibrium thickness of the layer 
becomes larger, giving the gas distribution a characteristic 
"flared" appearance. The exact form of this flaring depends 
on the amount of mass close to the plane of the Galaxy. 
Thus, by comparing the observed flaring to the predictions 
of the hydrostatic equilibrium arrangement of gas in the 
mass models developed in Section ^, we can see what degree 
of halo flattening is consistent with the observations. 



5.1 The observed flaring of the gas layer 

Before we can apply this technique, we need to summarize 
the observational data available on the flaring of the Galac- 
tic Hi layer. Merrifield (1992) calculated the thickness of 
the gas layer across a wide range of Galactic azimuths in 
his determination of the outer Galaxy rotation curve. As a 
check on the validity of that analysis, we have also drawn 
on the work of Diplas & Savage (1991) and Wouterloot et 
al. (1992), which derived the gas layer thicknesses across a 
more limited range of azimuths. For completeness, we have 
also included the data for the inner Galaxy as derived by 
Malhotra (1995). The resulting values for the full- width at 
half maximum (FWHM) of the density of gas perpendicular 
to the plane are given in the upper panel in Fig. ^. Note 
that, once again, the results depend on the adopted values 
of the Galactic constants - since the radii in the Galaxy of 
the various gas elements were derived from their line-of-sight 
velocities via equations (jjj) and (|^), the values of R in Fig. ^ 

depend on Ro and 0oQ 

It would appear from this figure that there are some dis- 
crepancies between the various measurements in the outer 
Galaxy. After some investigation, we established that these 
differences can be attributed to the effects of the beam sizes 
of the radio telescopes with which the observations were 
made. Such resolution effects will mean that the FWHM of 
the gas will tend to be overestimated. In the lower panel, 
we show what happens when the appropriate beam correc- 
tion is made to the data from Diplas & Savage (1991) and 
Wouterloot et al. (1992). No similar correction is required 
for the Merrifield (1992) analysis, as in that work the de- 
rived value of the FWHM was dominated by gas towards the 
Galactic anti-centre, which lies at relatively small distances 
from the Sun, and so the beam correction is small. Clearly, 



5 HALO FLATTENING FROM GAS LAYER 
FLARING 

We now turn to the technique developed by Oiling (1995) 
for measuring the shape of a dark halo from the observed 
thickness of a galaxy's gas layer. In essence the approach 
is similar to the stellar-kinematic method described above: 
the over-all mass distribution of the Milky Way is inferred 



t Merrifield's method for determining the thickness of the gas 
layer also exposes some of the shortcomings of more traditional 
methods. If the wrong values for Ro,@o and 0(R) are chosen, 
the inferred thickness of the gas layer (at a given R) will show a 
systematic variation with Galactocentric azimuth. It is possible 
to correct published data for this effect, but only if the assumed 
rotation curve, Galactic constants, as well as the thickness of the 
gas layer as a function of azimuth are specified. 
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Figure 4. The thickness of the Hi layer of the Milky Way. The widths in the inner Galaxy were taken from Malhotra (1995). For the 
outer Galaxy, we plot the widths from Diplas & Savage (1991; open triangles), Wouterloot et al. (1992; crosses), Merrifield (1992; open 
circles), and the average (thick full line & filled squares). The top panel represents the "raw" measurements, in bottom panel we present 
the beam-size corrected widths. 



this correction brings the various data sets into much closer 
agreement. We therefore adopt the mean curve shown in this 
panel for the following analysis; the error bars shown rep- 
resent the standard error obtained on averaging the various 
determinations. 



5.2 Sources of pressure support 

In order to compare the observed gas layer flaring to the pre- 
dictions for a gas layer in hydrostatic equilibrium, we must 
address the source of the pressure term in the hydrostatic 
equilibrium equation. The most obvious candidate for sup- 
porting the Hi layer comes from its turbulent motions. In 
the inner Galaxy, the Hi is observed to have a ve locity dis- 



)ersi on of o g — 9.2 kms 



independent of radius (Malhotra 



1995), and we assume that this value characterizes the tur- 



bulent motions of the gas throughout the Galaxy, providing 
a kinetic pressure term. 

Potentially, there may be other forces helping to sup- 
port the Galactic Hi layer: non-thermal pressure gradients 



associated with magnetic fields and cosmic rays may also 
provide a net force to resist the pull of gravity on the Galac- 
tic Hi. However, the analysis we are doing depends most on 
the properties of the Hi layer at large radii in the Milky Way, 
where star formation is almost non-existent, so energy input 
into cosmic rays and magnetic fields from stellar processes 
is likely to be unimportant. 

Our concentration on the properties of gas at large radii 
also eliminates another potential complexity. In the inner 
Galaxy, the interstellar medium comprises a complicated 
multi-phase mixture of molecular, atomic and ionized mate- 
rial. A full treatment of the hydrostatic equilibrium of such 
a medium is complicated, as gas can transform from one 
component to another, so all components would have to be 
considered when calculating hydrostatic equilibrium. For the 
purpose of this paper, it is fortunate that at the low pres- 
sures characteristic of the outer Galaxy it is not possible 
to maintain b oth the cold molecular phase and the warm 
atomic phase (Maloney 1993; Wolfire et al. 1995). Braun & 
Walterbos (1992) and Braun (1997,1998) have shown that 
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the cold phase disappears when the B-band surface bright- 
ness of a galaxy falls below the 25th magnitude per square 
arcsecond level, which occurs at ~ 1.5i?o in the Milky Way 
(Binney & Merrifield 1998, §10.1). Further, the ionized frac- 
tion of the ISM is expected to decrease with dist ance as it 
is close ly associated with sites of star formati on (Ferguson 



et al. 1 )96; Wang, Heckman & Lehneert 1997). Ultimately, 
the ionizing effects of the extragalactic background radia- 
tion field become significant, but onl y when the H i column 
density falls below about lAfp pc -2 ( Maloney 1993 ; Dove & 
Shull l! jj944| ; |011ing 1996c] ), which lies well beyond the radii 
we are considering here. 

For the current analysis, it therefore seems reasonable to 
treat the Milky Way's gas layer as a single isothermal com- 
ponent supported purely by its turbulent motions. More- 
over, this assumption has been made in the previous imple- 



mentat ions of the gas flaring meth od (Oiling 1996b 
quaert 



Combes & Viallefond 1997). A principal objective 



Bcc- 



of this paper is to test the validity of those analyses by com- 
paring results obtained by the flaring technique to those ob- 
tained by other methods. It is therefore important that we 
make the same assumption of a single isothermal component 
in the present study. 

5.3 Fitting to model gas layer flaring 

We are now in a position to compare observation and the- 
ory. The technique for calculating the gas layer thickness in 
different mass models has been described in detail by Oiling 
(1995). In brief, for each model, at each radius R, one inte- 
grates the hydrostatic equilibrium equation, 



9$ 

dz 



1 dp g 



Pa 



dz 



(T) 



to obtain the gas density distribution, p g (R, z). The FWHM 
of this model gas distribution can then be compared directly 
with the observations. 

The results of such calculations are illustrated by the 
examples shown in Fig. [|. The basic trends in this analysis 
are clearly demonstrated by these examples. The decrease 
in total density with radius leads to a dramatic flaring in 
the model-predicted gas layer thickness, just as is seen in 
the observations. For a flatter model dark halo, the mass is 
more concentrated toward the plane of the Galaxy, squeezing 
the Hi layer into a thinner distribution. 

Once again, the results depend quite sensitively on the 
choice of Galactic constants, since these values affect both 
the gas distribution as inferred from observations and the 
acceptable mass models as inferred from the rotation curve. 
As is apparent from Fig. none of the Ro = 8.5 kpc models 
fits the observations. In fact, to match the observed layer 
width one would require a substantially prolate dark matter 
halo with q ~ 1.5, which none of the current dark matter 
scenarios predict. For Rq = 7.1 kpc, on the other hand, a 
very good fit is obtained for models with a halo flatness of 
q ~ 0.7. Such models even reproduce the observed inflection 
in the variation of the gas layer width with radius at R ~ 
10 kpc. 

For each plausible set of Galactic constants, we can 
carry out a similar analysis to that in Section ^ to see what 
range of values of q are consistent with the observed gas 
layer flaring. Because this technique relies on data from large 




5 10 15 

Galactocentric radius (kpc) 



Figure 5. The observed flaring of the Galaxy's gas layer (trian- 
gles) compared with the flaring predicted for a number of mass 
models that reproduce the observed rotation curve and value of 
E*. The models in the top panel have been calculated assum- 
ing Ro = 8.5 kpc and a disk mass-to-light ratio T d k— band = 
0.60; the bottom panel shows models with Rq = 7.1 kpc and 
Yd K-band = 0.41. The models shown here have a disk scale- 
length of hd = 2.5 kpc. The halo flattening of the model and the 
reduced x 2 value of the fit is shown for each case. 



radii in the Galaxy, where the dark matter halo is the dom- 
inant source of mass, the flaring predicted by the models 
depends very little on the properties of the disk and bulge. 
Unlike the stellar kinematic analysis, therefore, one cannot 
trade off the mass in the disk against the mass near the 
plane from a more-flattened dark halo. This difference is il- 
lustrated in Fig. |§|, which shows the way that the value of 
q inferred from the gas layer flaring depends on the prop- 
erties of the stellar disk (as parameterized by the model's 
column density of stars in the Solar neighbourhood). As for 
the stellar-kinematic analysis, there is a well-defined correla- 
tion between q and £» , but, for the reasons described above, 
the trend for the current method is very much weaker, and, 
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within the observationally-allowed range for £*, q is tightly 
constrained to 0.73 ± 0.03. 

Although this constraint is remarkably good, it should 
be borne in mind that it is still dependent on the adopted 
values for the Galactic constants. As we have seen above, 
larger values for Ro lead to rounder, or even prolate esti- 
mates for halo shape, so we will not obtain an unequivocal 
measure for q from this method until the Galactic constants 
are measured more accurately 



6 COMBINING THE TECHNIQUES 

The different slopes of the two relations in Fig. || raises an 
interesting possibility. Clearly, for a consistent picture, one 
must use a single mass model to reproduce both the stellar- 
kinematic constraint on the mass in the solar neighborhood 
and the observed flaring of the gas layer. Thus, although 
there are whole linear loci in this figure of models with dif- 
ferent values of q and E* that satisfy each of these con- 
straints individually, there is only the single point of inter- 
section between these two lines where the model fits both 
the stellar-kinematic constraint and the observed flaring of 
the gas layer. Hence, for given values of Ro and 0o, one 
predicts unique values for q and £,. 

We have therefore repeated the analysis summarized in 
Fig. ^ spanning the full range of plausible values for _Ro and 
0o, and calculated the mutually-consistent estimates for E* 
and q for each case. In order to reduce the computational 
complexity of this large set of calculations to manageable 
proportions, we made use of Olling's (1995) fitting formula, 
which showed that, if self-gravity is negligible, one can ap- 
proximate the model-predicted thickness of the gas layer by 
the relation 

FWHM(i?) « J 13 d 5q ° B y/RjTlP, (8) 
y 1.4 + q Vh^ 

where Vh,oo is the circular rotation speed of the dark halo 
component at large radii. Allowing for the self-gravity of 
the gas layer, one obtains a similar formula. Comparing the 
values derived from this formula to sample results obtained 
by the full integration process described above, we found 
that the approximation based on equation (^) matches the 
detailed integration for R > 1.75i?o- We therefore only used 
the data from beyond this radius in the following analysis, 
in which the model gas layer thickness was estimated using 
the approximate formula. 

The values obtained for q and E* for each possible pair 
of Galactic constants are presented in Fig. ^| Thus, for ex- 
ample, for values of Ro = 7.1 kpc and 0o = 185 kms -1 , one 
finds E* = 35Mq pc -2 and q = 0.7, corresponding to the 
intercept that we previously calculated in Fig. |j| 

Figure |^ places some interesting limits on the properties 
of the Milky Way. For example, if we maintain our prejudice 
that the halo should be oblate (q < 1), then, unless we adopt 
a particularly extreme value for _Ro, we find that 0o must be 
less than ~ 190 km s _1 . If we also adopt Kuijken & Gilmore's 
(1989) measurement of E* = 35 ± 5Mq pc~ 2 , we find that 
only models within the heavily-shaded region of Fig. [] are 
acceptable, placing an upper limit on _Ro of ~ 7.6 kpc. Con- 
versely, if one forces the Galactic constants to the IAU stan- 
dard values of Ro = 8.5 kpc and 0o = 220kms _1 (Kerr 



& Lynden-Bell 1987), one finds barely-credible values of 
£, ~ 6OM pc -2 and q ~ 1.5. 



7 CONCLUSIONS 

The prime objective of this paper has been to check the 
validity of techniques for measuring halo flattening by ask- 
ing whether two different techniques return consistent values 
when applied to the Milky Way. As we have seen in the last 
section, the answer is a qualified "yes." The qualification is 
that consistency with the measured stellar column density 
requires values of the Galactic constants that differ from 
those conventionally adopted. However, as was discussed in 
Section |^, the true values of these constants remain elusive, 
with estimates spanning the ranges 7 kpc < Ro < 8.5 kpc 
and 185kms _1 < 0o < 235kms _1 . With such gross uncer- 
tainties, it quite straightforward to pick values that produce 
an entirely self-consistent picture. To underline this point, 
we have included on Fig. |^ the results of Oiling & Merrifield's 
(1997) estimates for Ro and ©o derived from an analysis of 
the Oort constants. If that analysis is valid, then we have a 
consistent model for the Milky Way in which the dark halo 
has a flattening of q ~ 0.8. 

Ultimately, this analysis will allow us to come to one of 
two conclusions: 

(i) If future studies confirm low values for Ro and 0o sim- 
ilar to those derived by Oiling & Merrifield (1997), then the 
consistency of the two analyses for calculating q imply that 
the gas layer flaring technique is valid, adding confidence to 
the previous determinations by this method. 

(ii) Conversely, if we learn in future that Ro and ©o are 
closer to the more conventional larger values, then the im- 
plied values of E* are so far from the observed estimates that 
one has to suspect that at least one of the techniques for 
measuring q is compromised. In this case, one would have to 
look more closely at some of the assumptions that went into 
the analysis. For example, perhaps the non-thermal pres- 
sure forces from cosmic rays and magnetic fields have a role 
to play even at large radii in galactic disks. Alternatively, 
perhaps the Hi layer is not close enough to equilibrium for 
the hydrostatic analysis to be valid. Finally, our assumption 
of azimuthal symmetry may be invalid. Strong departures 
from axisymmetry could mean that our determination of the 
thickness and column density of the gas is compromised, and 
that the locally determined values of Ei.i and E* may not 
be representative for the Galactocentric radius of the Sun. 

Assuming for the moment that the analysis is valid, we have 
another datum to add to Fig. |l| Since the two previous flar- 
ing analyses returned systematically rather small values of 
q ~ 0.3, it is reassuring that the Milky Way seems to indi- 
cate a larger value of ~ 0.8 - it appears that the low values 
are simply a coincidence arising from the very small number 
statistics. This larger value is inconsistent with the very flat 
halos that are predicted by models in which the dark matter 
consists of either decaying neutrinos (Sciama 1990) or cold 
molecular hydrogen (Pfenniger, Combes & Martinet 1994). 

With the addition of the Milky Way to the data pre- 
sented in Fig. 0, the only technique that stands out as giv- 
ing systematically different values for q is the bending mode 
analysis of warped gas layers. In this regard, it is notable 
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Figure 6. Contours of the stellar column density in the Solar neighbourhood long dashed lines) and halo flattening (q; dotted lines) 
as a function of the adopted values for the Galactic constants, Ro and ©o- The heavy full line and the heavy dashed line corresponds to 
Kuijkcn & Gilmore's (1989) determination of E«, and the ±1 — <r values. The heavily-shaded region corresponds to parts of parameter 
space consistent with these values of S* that produce an oblate halo. The lightly-shaded area gives the corresponding region if we adopt 
Gould, Bahcall and Flynn's (1997) values for S t . The cross shows the determination of the Galactic constants based on an analysis of 
the Oort constants (Oiling & Merrifield 1998). 



that simulations cast some doubt on the validity of such 
analyses. The method assumes that warps are manifesta- 
tions of persistent bending modes that occur when the flat- 
tening of a disk and the surrounding dark halo are mis- 
aligned. However, the simulations show that gravitational 
interactions between a misaligned disk and halo rapidly 
bring the two back into alignment, effectively suppressing 
this mechanism (e.g. Dubinski & Kuijken 1995, Nelson & 
Tremaine 1995, Binney, Jiang & Dutta 1998). The number 
of measurements is still rather small, but it is at least inter- 
esting that if one neglects the warped gas layer results, the 
remaining data appear entirely consistent with the dotted 
line in Fig. Ill which shows Dubinski's (1994) prediction for 



the distribution of halo shapes that will be produced in a 
cold dark matter cosmology. 
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